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"second-class" hadronic matrix elements must be determined in order to exact such tests. 
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I. INTRODUCTION 

Inspired by the pioneering global fits of the CKM ma- 
trix developed by the CKMfitter Group [1, 2] and the 
UTFit Collaboration [3-7] for the interpretation of flavor- 
physics results from the B factories and the Tevatron, 
wc outline the prospects for the elucidation of physics 
beyond the standard model (BSM) via a global fit of 
neutron ^ decay observables, including the lifetime and 
the energy-dependence of the angular correlation coeffi- 
cients. Our global fit, which we term nFitter, employs 
a maximum likelihood statistical framework which ac- 
counts for both experimental and theoretical uncertain- 
ties, with the latter arising primarily from the poorly 
known weak hadronic second-class currents. 

The elucidation of a 'V — A" law [8, 9] in the media- 
tion of low-energy weak interactions played a crucial role 
in the rise of the standard model (SM) [10]. A variety 
of wcll-motivat(>d SM extensions speak to the possibility 
of new dynamics at the Fermi scale, which we identify 
via V = {2\/2Gf)~^^'^ ~ 174 GeV, and, concomitantly, 
to tangible departures from the V — A law in low-energy 
weak processes, be it through scalar or tensor interac- 
tions, or right-handed currents. Alternatively, e.g., new, 
light degrees of freedom could appear and yield viola- 
tions of the V — A law — if probed at sufficient exper- 
imental resolution. At the same time, an ongoing vig- 
orous experimental program for precision measurements 
of neutron /3 decay observables with cold and ultracold 
neutrons [11 17] exists with an overarching goal of real- 
izing bettered assessments of the limits of the SM. The 
experimental effort focuses on measurements of two gen- 
eral types of observables: angular correlation coefficients, 
which parametrize the angular correlations between the 
momenta of the various decay products and/or the spin 
of the initial state neutron or electron in the differential 
decay rate, and the neutron lifetime, r. Of the future ex- 
perimental plans reviewed in Refs. [11-19], the claim of 
ultimate precision rests with the PERC experiment [20], 
for a sensitivity of up to 10~* precision. 

The angular correlation coefficients a, b, A, B, and 
D [21], of which only a, A, and B are appreciably 



nonzero [22], describe the distribution in electron and 
neutrino directions and electron energy in neutron de- 
cay. Since the nucleon mass is markedly larger than the 
neutron-proton mass difference, a recoil expansion of the 
differential decay rate of the neutron in terms of the ra- 
tio of various small energy scales to the nucleon mass is 
extremely efficient. The a, A, and B correlation coeffi- 
cients in this order, which we denote with a "0" subscript, 
arc functions only of A = qa/qv > in the SM, where 
qa and gv are the weak axial- vector and vector coupling 
constants; this is tantamount to the V — A law. Thus, 
extractions of gq, Aq, and Bq from the measured param- 
eters a, A, and B determine A in the SM, which, by it- 
self, is a fundamental parameter of the weak interaction. 
Using the measured values for t and A determines gv 
and gA independently. The former, once radiative correc- 
tions are applied [23, 24], yields the Cabibbo-Kobayaslii- 
Maskawa (CKM) matrix element Vud- A neutron-based 
value for Vud [22] is not yet competitive with the result 
extracted from measurements of the ft values in super- 
allowed 0+ — > 0+ nuclear /3 decays [25, 26]. Such efforts 
are nevertheless well-motivated in that such is not sub- 
ject to the nuclear structure corrections which must be 
applied in the analysis of 0+ — > 0+ transitions. We set 
such prospects aside and focus on the possibility of dis- 
covering dynamics beyond the SM through an assessed 
quantitative violation of the V — A law. 

We consider neutron beta decay observables exclu- 
sively, though the framework — and fits — we employ can 
readily be enlarged to encompass all low-energy semilep- 
tonic processes involving first-generation quarks. This is 
possible because we can employ, rather generally, a sin- 
gle, quark-level effective Lagrangian, at leading power in 
v/A [27, 28], for all such processes [29, 30], where we re- 
fer to Ref. [18] for a review. To realize this we need only 
assume that new physics appears at an energy scale A 
in excess of the scale v, crudely that of the W and Z 
masses. This allows the construction of an effective La- 
grangian in terms of operators of mass dimension A*^"^ 
with d > 4 [31], where the nonobservation of non-SM in- 
variant operators in flavor physics [2, 7, 32], allows us to 
impose SM electroweak gauge invariance in its construe- 
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tion. At leading power in u/A, there are then precisely 
ten dimension-six operators describing the semileptonic 
decay of the d quark. The coefficients of these opera- 
tors can be determined, or at least constrained, through 
a global fit of (3 decay observables — in neutrons and 
nuclei — and of meson decay observables as well. This 
quark-level description, upon matching to a nucleon-level 
effective theory [33] , yields a one-to-one map to the terms 
of the effective Hamiltonian constructed by Lee and Yang 
[34] , admitted by Lorentz invariance and the possibility of 
parity violation [34] . The latter framework is employed in 
the usual analysis of the angular correlations in beta de- 
cay [21]. The construction of the underlying quark-level 
effective theory allows the inclusion of meson observables 
on the same footing. 

Global fits of the Lee- Yang coefficients have been made 
to beta decay observables, particularly of the phase- 
space-intcgrated angular correlation coefficients [12, 35]. 
As we develop explicitly, the advantage of the current 
approach is that we can take theory errors into direct 
account in the optimization procedure. We focus on the 
simultaneous fit of the angular correlation coefficients a 
and A in neutron decay and of the neutron lifetime to 
limit the appearance of tensor and scalar interactions 
from physics BSM, in part because the theory errors in 
this case are spanned by the ill-known nucleon matrix 
elements of particular local operators, which presumably 
and eventually can be computed in QCD using the tech- 
niques of lattice gauge theory. 

The empirically determined correlation coefficients and 
neutron lifetime, fixed to some precision, can be used to 
set limits on new dynamics. In this regard it is natural to 
focus on the Fierz interference term b because it is linear 
in new physics couplings. That is, it is possible to dis- 
cover evidence of scalar or tensor interactions, through a 
measurement of b in excess of the ^ 10""^ level expected 
in the SM from recoil-order effects [36]. To date there 
have been no published results for b in neutron /3 decay, 
although several efforts are underway. Searches for scalar 
and tensor interactions have also been pursued via mea- 
surements of the decay electron's transverse polarizations 
with respect to the neutron spin [37, 38] , though no new 
experiments of such ilk arc currently planned. Informa- 
tion on b can be gleaned in different ways. Its presence 
modifies the shape of the electron energy spectrum in the 
differential decay rate, and direct searches probe that. It 
is also subsumed in measurements of the energy depen- 
dence of the a, A, and B correlation coefficients [33, 39], 
so that we wish to consider the additional experimental 
observables offered by (electron) energy dependence with 
some care. The latter pathways offer "indirect" access 
to b and require greater statistical control than spectral 
shape measurements for fixed sensitivity to scalar and 
tensor interactions. However, they arc also less sensitive 
to systematic errors, both experimental and theoretical. 
In this first paper we focus on access to scalar and tensor 
interactions through measurements of a, A, and r, as we 
now explain. 



The V ~ A description of neutron (3 decay includes 
contributions from six possible weak hadronic currents 
which give rise, at recoil order, to energy-dependent ex- 
pressions for the angular correlation coefficients. How- 
ever, the most recent extractions of A from A [40-44], 
which yield, at present, the most precise determination 
of A, have assumed the validity of the Conserved- Vector- 
Current (CVC) hypothesis (i.e., the CVC value for the 
weak magnetism coiipling) and neglected second-class 
currents. This is certainly reasonable at the current level 
of experimental precision, although measurements of A at 
the ^ 0.5% level of precision are, in principle, sensitive to 
the weak magnetism coupling, which contributes to the 
asymmetry at the ^ 1.5% level. With the anticipated 
increased sensitivity to a, A, and B in next-generation 
experiments, we note that precise measurements of the 
energy dependence of the a and A correlations offer the 
possibility to test CVC and to search for second-class cur- 
rents independently [45]. In nFitter we fit the energy 
dependence of the angular correlation coefficients directly 
and thus need not assume either the validity of the CVC 
hypothesis or the absence of second-class currents, yield- 
ing a framework for a robust test of the validity of the 
V — A description of neutron j3 decay. 

The outline of the remainder of this paper is as follows. 
First, we review the formalism for the analysis of neutron 
j3 decay observables in Section II. We then briefly review 
in Section III the experimental status of, or limits on, the 
various SM parameters relevant to our global fit, such as 
second-class currents. We then describe the maximum 
likelihood approach to our global fit in Section IV, where 
we discuss the construction of our likelihood function, as 
well as the inclusion of experimental and theoretical un- 
certainties. We illustrate the prospects of our global fit 
with a few numerical examples employing the frequentist, 
or specifically i?Fit [1, 2], statistical procedure, reserving 
the use of alternative statistical procedures for later work. 
We then show examples of the results from nFitter fits 
to Monte Carlo-generated pseudodata in Section V un- 
der various scenarios, and we quantify via these exam- 
ples the statistical impact that future improvements in 
the precision of neutron 13 dcicay observables should have 
on the assessment of the validity of the SM. We also dis- 
cuss the extent to which theoretical uncertainties in the 
presently poorly known second class contributions limit 
such assessments, and we can extract from such stud- 
ies the precision to which they should be established to 
obviate that impact. Finally, we conclude with a brief 
summary in Section VI. 



II. FORMALISM FOR NEUTRON /3-DECAY 
OBSERVABLES 

If we suppose the low-energy, effective weak interac- 
tion, employing explicit n and p degrees of freedom, is 

mediated by the ten dimension-six operators enumerated 
by Lee and Yang [34] , then the differential decay rate for 
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neutron /3 decay takes the form [21] 
dT 1 
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where we refer to Ref. [21] for the cxphcit form of ^ and 
the correlation coefficients in terms of the parameters 
of the Lee- Yang Hamiltonian [34], noting Ref. [18] for 
a discussion of the connection to modern conventions. 
We use Ef, {E„) and p^ {pv) to denote, respectively, the 
electron's (antincutrino's) total energy and momentum, 
where Eq is the electron endpoint energy, and ((?„) is the 
neutron polarization. 

The Coulomb corrections to Eq. (1) are also known 
[46] and modify the expression most notably in terms of 
a multiplicative Fermi function F{Z, E^) [47]. The phase- 
space integrated Fermi function and corrections to it have 
been studied in great detail [48, 49]; we omit it, as well 
as the outer radiative correction [50], in the generation 
of the Monte Carlo pseudodata for our decay correlation 
studies as we are interested in a{Ee) and A{Ee), which 
are accessed through asymmetry measurements for which 



such effects only lead to a slight modification of the rel- 
ative statistics (via the spectral shape). 

The D term is a naively time-reversal-odd observable: 
a value for D in excess of the ~ 10~^ level attributed 
to SM final-state interaction effects [51, 52] would reveal 
the existence of new CP-violating interactions at the La- 
grangian level (assuming CPT holds). The current level 
of experimental precision places stringent constraints on 
any such new effects [53, 54]. 

In what follows we report expressions for the correla- 
tion coefficients which include the trcc-lcvcl new physics 
of the Lee- Yang Hamiltonian and the contributions of the 
usual V — A terms through recoil order. In realizing this 
the strong interaction plays an essential role: the ma- 
trix elements of the vector V and axial- vector A currents 
are described by six distinct form factors. We find it 
immensely useful to note the quark-level effective theory 
which underlies the Lee- Yang couplings [29, 30, 33]; such 
makes the separation of the QCD physics which underlies 
the hadronic matrix element calculation from the nomi- 
nally higher-energy physics encoded in the effective low- 
energy constants clear. As per Refs. [30, 33] we map the 
Lee- Yang effective couplings Ci,C[ with i G {V, A, S,T} 

to cP = {Gf / V^)VudCl'^ and note the hadronic matrix 
elements needed in (3 decay are [55] 
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Unip) , 



{p{p')\u-f''d\n{p)) = up{p') 

{p{p')\u^^l5d\n{p)) = Up{p') 

{p{p')\ud\n{p)) = Up{p')gs{q^)un{p) , 
W)\ucr^ud\n{p)) ^ upip') [gT{q')<y'"' + g^T\<i^)i<i''r - q'ln 

+9P{q^WP'' - q^Pn+g^Til^KrH" - YHn] Unip) , 



(2) 

(3) 

(4) 

(5) 



where q = p' — p denotes the momentum transfer, P = 

p'+p, and M is the neutron mass. In neutron /3 decay, the 
g^-dependent terms are of next-to-next-to-leading order 
(NNLO) in the recoil expansion, noting /i(0) and gi{0) 
appear in leading order (LO), and hence are of negli- 
gible practical relevance. Consequently, we replace, as 
usual, the form factors with their values at zero momen- 
tum transfer. We note /i(0) = gv is the vector coupling 
constant given by (/y = 1 under CVC; /2(0) = /2 is the 
weak magnetism coupling constant given by {Kp — k„)/2 
under CVC, noting Kp(„) is the anomalous magnetic mo- 
ment of the proton (neutron); /3(0) = /s is the induced 
scalar coupling constant; gi{0) — gA is the axial vec- 
tor coupling constant; 52(0) = 52 is the induced tensor 
coupling constant; and 53(0) = gp is the induced pseu- 
doscalar coupling constant. The CVC predictions have 
SM corrections in NNLO. The contributions of fi, f2, gi, 



and 33 to the hadronic current are termed first-class cur- 
rents, whereas those of /a and 92 are termed second-class 
currents, due to their transformation properties under 
G-parity [55]. The latter quantities, /3 and g2, vanish 
in the SM up to quark mass effects which break flavor 
symmetry; we discuss their estimated size in Sec. 111. 

Of particular interest to us are the scalar and tensor 
interactions, as establishing their existence at current ex- 
perimental limits would signify the presence of physics 
BSM. The matching of the quark-level to nucleon-level 
effective theories at LO in the recoil expansion yields: 

Cs = gsif^s + es) , 

C's = gsies - es) , 

Ct = 4gT(eT + er) , 

C't = ^grier - ?t) , (6) 
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where the e coefRcients are the low-energy constants of 
the quark- level effective theory of Refs. [30, 33]. We have 

(i) 

neglected the matrix elements g)p' with i e 1, 2, 3 in re- 
alizing this expression and thus, for consistency, shall 
neglect the scalar and tensor contribution to recoil order 
terms in all that follows. Bhattacharya et al. [33] have 
employed a iJfit scheme to determine the impact of im- 



proved lattice estimates of gs and gx on the limits on the 
quark-level low-energy coeffcients for given experimental 
sensitivities to Cg^j,. 

In unpolarizcd neutron /3 decay, the unpolarized dif- 
ferential distribution relevant for a measurement of a, 
neglecting terms beyond next-to-leading order in the re- 
coil expansion but accounting for all six possible form 
factors, is of the form [56] 
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( 4A(1 + Kp- Kn) 
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A52 



l-h3A2 



(7) 



where Ogv is the electron-antineutrino opening angle and /3 = \pe\/E(.. The structure of this expression serves as a 

de facto definition of a = ai + a2f3 cosOeu and &bsm in recoil order. It follows that of Rcf. [21] if recoil terms are 
neglected and is that of Ref. [58] if &bsm = 0. Note that in writing the recoil contributions we have neglected terms 
of 0{esgs, ^TgT)R- Moreover, 



ai = ao + 



1 



(1 + 3A2)2 



4A(1 + A + A2 + A^ + 2/2 + 2/2A2)i? + (1 -I- 2A - - A^ + 4/2A - Af^X^) 



- [8A(1 + 2/2 + A2 + 2/2A2) + 3(1 + 3A2)2] Rx 
+ [2(A - A3)52 + 2(A2 - l)h]^ + 8A(1 + \^)g2R 



Rx 



3(A2 - 1) ^ 



(1 + 3A2)' 

with A = gA/dv > in the SM, and the kinematic factors e, R, and x are defined according to 



/ m-e \ 2 

Km) ' 



R = 



Eq 
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X = 



Ee 
En 



(8) 
(9) 

(10) 



The computations of Ref. [57] have been repeated in de- 
riving these forms, and the results arc consistent up to 
the /a terms [45]. It is also consistent with Ref. [36], as 
well as with Ref. [58], noting f^ = g2 = 0m the latter. 
These comparisons are all within the context of V — A 
theory. 

Wc use R itself, noting R « 1.37 x 10^^, to characterize 
the efficacy of the recoil expansion. Both SM and BSM 
couplings appear in S, aq, and &bsm, namely [21] 



S = 1 + 3A2 + {gsesf + SiAgrerf 



(11) 



„ (l-A2)-(gg65)2 + (4ffTeT)^ 



&BSM 



I + 3A2) + {gsesr + K^gTeT^ 

2(.gge5) - 6A(4gT£T) 
(l + 3A2) + (gses)2 + 3(4ffTeT)2 



(13) 



where we employ Eq. (6). 

Our recoil-order expression for the term proportional 
to e/Rx oc nie/Ef. appearing within the first set of 
square brackets in the differential distribution of Eq. (7) 
is equivalent to the term labeled "5sm" employed in Refs. 
[33, 59]. However, it should be noted that the second- 
class currents /s and g2 yield an additional me/ Eg term 
which is proportional to (/a — Ag2)- Simply for the sake 
of notation, we label this term "6scc"i where we then 
have, in summary. 



me 1 -I- 2A A2 -I- 2X{Kp - k„) 



&SCC = 2 



M l-h3A2 



(14) 



In polarized (3 decay, the differential distribution rele- 
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vant to a measurement of A is of the form 

oc M^R^l3x\l - xf-, 

dE^(M^ ^ ^ ' (l + e-2i?a;)' 



1 + &BSM ^ + A/9 cos 6*6 



(15) 



where is the angle between the momentum of the elec- 
tron and the polarization of the neutron. Here, too, the 
structure of this expression follows that of Ref. [58] and 
Ref. [21] in suitable limits and serves as a definition of A 



in recoil order; note that wc neglect recoil contributions 
to 6bsm " so that "6bsm" is the same quantity hero and 
in Eq. (7). The complete expression for /t(.T) can be found 
in Ref. [56] and is in agreement with Ref. [60]. To LO in 
the recoil expansion, h{x) is of the form 



(16) 



Working to LO in the S and T terms and to NLO in the 
V — A terms, 



h{x) 
(1 + e - 2Rx)^ 



(1 + 3A2) 



1 I qff^ I p^/^ 4A(l + /Cp-K„) \ / A^ + A + A(Kp-/t„) 
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1 + 3A2 / Rx 
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Rx\l + 3A2 



+ (A^ -1)[Rx 
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+ (gsesf + H^^grerf 



(17) 



and A is of the form 



A = Ao 



(1 + 3A2)2 I Rx 



[4A2(1 - A)(l + a + 2/2) + 4A(1 - A)(A<?2 - /s)] 



with 



+ R 
+ Rx 

An = 



;[1 + A + 2(/2 + <?2)](3A2 + 2A-1) 



^(1 + A + 2/2)(l - 5A - 9A2 - 3A3) + ^^2(1 + A + SA^ + SA^) 



2A(1 - A) + 2{AgTeTf + 2{gses){^9TeT) 



(18) 



(19) 



(1 + 3A2) + {gsesY + 3(4gTeT)2 
Our expressions in the context of V—A theory agree with those of Ref. [36] — and with those of Ref. [58] if f3=g2 = 0. 



As a final topic we revisit the computation of the neu- 
tron lifetime and focus particularly on the role of recoil- 
order corrections. In the current state of the art [23, 24], 
Vud, 2, and r are related by 



r = 



4908.7(1.9) s 



(20) 



where S = 1+3A2 in the absence of new physics. Employ- 
ing Vud = 0.97425 and A = 1.2701 [22] yields a lifetime 
of 885.6 s. The numerical value reported in Eq. (20), as 
per Ref. [24] , incorporates an improved treatment of elec- 
troweak radiative effects, including certain 0{a^) contri- 
butions [23, 24]. Note that the calculation embeds a value 
of gA = 1-27 in matching the short- and long-distance 
radiative corrections [23, 24]. The numerical value also 
includes the phase space factor / which incorporates the 
Fermi function and various recoil-order terms [48]; many 
small terms involving hadronic couplings other than gy 
and gA can enter in recoil order. Reference [48] analyzes 



neutron /3 decay to 0.001% in precision and finds the lat- 
ter corrections negligible save for the possibility of that 
from 52- In that work [48] contributions proportional 
to /a are assumed to be strictly zero from CVC, though 
such can also be engendered by SM isospin violation. We 
can evaluate the various small contributions to the total 
decay rate by integrating the terms of Eq. (17) over the 
allowed phase space of Eq. (15). We denote a contribu- 
tion relative to that from S by Cg^g^, where gi and gj 
are the couplings it contains. Defining Wq = E^'^/rUe, 
W = Ee/nie, pw = VW^ — 1, and finally 



Im = dWpwWiWo-WfW"^ , (21) 



the small contributions to the decay rate, relative to that 
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mediated in LO by S, are 

Xg2 m. 



n — ^ 



^gvh 
CgAh 



2/3 me /_! 

S M 7o ' 
A(«;p - K„) me 4Ji - 2Wo/o - 2/_i 



S M V /o 
A me Uh-2WolQ-2I_i 
S M V io 
A^me / 2W0J0 + J-1 
S M I /o 



(22) 



where wc note the appendix of Ref. [48] for a useful tab- 
ulation of the integrals Im- In our current study, in 
which we assume that the entire range of possible elec- 
tron energies is experimentally accessible, both Cg^f^ and 
Cgyg^ vanish up to contributions nominally of 0{aR) ~ 
1.0 X 10~5 and ^ I.9 x lO"'^ in size, both of which 

are negligible at 0.001% precision. Using the masses re- 
ported in Ref. [22], specifically M = 939.565379 MeV, 
M' = 938.272046 MeV, and rUe = 0.510998928 MeV, the 
remaining terms evaluate to 



Cg^g, = -6.21 X 10-3^ 

=7.12x 10-*4, 



a 



gAgA 



\2 

-3.11 X 10"^— . 



Using A = 1.2701 [22] we find C, 



gAgA 



-8.58 X 10- 



If A changes within ±0.10, we note that Cg^g^ changes 
negligibly at the precision to which we work, so that we 
can regard Cgj^g^ as a fixed constant in the optimizations 
to follow. This particular contribution should already be 
embedded in the numerical constant of Eq. (20); however, 
the terms involving second-class currents have not been. 
To include such small corrections in the lifetime we need 
only replace S with S(l + Cg.g.), so that to retain 52, 
e.g., we write 



4908.7 



r = 



F„rf|2(S-(6.21xl0-3)52A) 



(24) 



It is worth noting that finite experimental acceptance 
plays an important role in the assessment of the recoil 
corrections to the lifetime. For example, if the accessi- 
ble electron kinetic energy were limited to the interval 
[100, 700] keV from the allowed range of [0, E^^"" -m^K, 
781.5] keV, then and Cg^gj^ would no longer vanish 

in 0(if!). The integrals are still analytically soluble and 
evaluate to C^^/, = 7.47x 10"* and Cg^g^ = 2.02x10"'*, 
where we use /2 = (Kp — Kn)/2 = 1.8529450 [22] as 
well. Taken together, they yield a contribution some 100 
times larger than our earlier assessment, which was set by 
0{aR). Including these effects in the manner of Eq. (24), 
they reduce the determined neutron lifetime for fixed A 



by 0.8 s to yield 884.8 s. Such considerations differenti- 
ate neutron lifetime experiments which (i) count surviv- 
ing neutrons from those which (ii) count decay products. 
The existing tension between the latter, "in beam" ex- 
periments and the former, "bottle" experiments [17, 22] 
— though there is also tension between the results of the 
most precise bottle experiments [17, 22] — make the ob- 
servation intriguing. The most precise in-beam neutron 
lifetime experiment [61, 62] counts decay protons, rather 
than electrons, so that our numerical analysis is not di- 
rectly relevant. However, analogous recoil contributions 
impact the proton energy spectrum as well [63] and simi- 
larly certain recoil terms also become of 0{R) rather than 
of 0{aR, i?^), once the finite experimental acceptance is 
included. This issue deserves detailed study. 



III. SURVEY OF THEORETICAL 
UNCERTAINTIES 



In the SM the corrections to the predictions of the 
CVC hypothesis, gy = 1, ./2 = (^p — Kri)/2, and /s = 
are parametrically known to be of 0{m,d — niuf [64] for 
gv and of 0{md — niu) for /2 and /a. The coupling 52 
does not vanish under CVC, but rather from its G-parity 
properties; it, too, is nominally nonzero at 0{md — m.„). 

The CVC hypothesis is assumed (i.e., gv = 1) in ex- 
tracting values for Vud from measured ft values in super- 
allowed 0+ 0+ nuclear /3 decay, and the universality 
(23) of 9v in these decays has been tested to 1.3 x 10""' at 
68% confidence level (CL) with a concomitant constraint 
of me/s/Afgy = -(0.0011 ± 0.0013) [25], implying h is 
constrained to only 0(1). Direct computation reveals the 
deviation of gv from unity to be smaller still [65]. De- 
cay correlation measurements invariably confiate tests of 
the CVC prediction for the weak magnetism form factor 
[66] with those which would limit second-class currents; 
currently the CVC value of /2 is tested to the level of 
some 6% [12, 67] at 68% CL. We note, however, there 
has not, to date, been a published measurement of /2 in 
neutron /3 decay, such as, e.g., could be extracted from 
the linear energy dependence of o, or ^ if 32 = is as- 
sumed. The second-class couplings and 92 have also 
not been probed experimentally in neutron P decay. The 
comparison of ft values in mirror transitions can also test 
for second-class currents; in that context the weak mag- 
netism form factor does not enter but an isospin-breaking 
additive correction from the axial form factors can. Such 
experiments give the strongest empirical constraints on 
second-class currents [68, 69] although additional theo- 
retical uncertainties enter. A survey of nuclear f3 decay 
data gives the limit l5'2//2l < 0.1 at 90% CL, yielding 
I52I < 0.2 at 90% CL [70]. 

Some theoretical studies of 92 exist, particularly in the 
case of strangeness-changing transitions. A bag model 
estimate gives g2/gv ~ 0.3 [71] in |A5| = 1 semilep- 
tonic transitions. More recently, non-zero second-class 
currents have been observed in quenched lattice QCD 
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calculations of form factors which appear in the hyperon 
semileptonic decay — )• yielding fz/gv — 0.14(9) 

and gild A = 0.68(18) [72]. Turning to the nucleon sec- 
tor, we expect these estimates to be suppressed, crudely, 
by md/ms ~ 0.1. This makes them nearly compatible in 
scale with the value for g2 determined using QCD sum 
rule techniques, g2/gA = -0.0152 ± 0.0053 [73]. 

The same lattice study has explored SU(3) break- 
ing in the /2 coupling as well, finding [/2//i]30^s+ ~ 
1.16(11) X [/2//i]n->-j)i a result somewhat different from 
the predictions of common models of SU(3) breaking [72]. 
Applying a scaling factor of md/ms ^ 0.1 we would sup- 
pose that CVC-breaking in /2 is no larger than a few 
percent. It is notable that the experimental limits on /2, 
/a, and g2 arc all rather lax with respect to theoretical 
expectations of CVC breaking and SCC from SM physics. 

In what follows we explore the impact of a non-zero 
52, as well as of values of /s and /2 which are not fixed 
precisely by the CVC prediction. These form factors all 
appear in recoil order; the usual first assumption of un- 
correlated errors suggests that the impact of these form 
factors on the potential discovery of BSM physics ought 
be modest. This turns out to be not so because the fit 
parameters, rather, are highly correlated, as we shall see. 



IV. MAXIMUM LIKELIHOOD ANALYSIS 
A. Construction of the Likelihood Function 

Having reviewed the formalism for neutron (5 decay, 
the starting point for our maximum likelihood analysis 
of neutron j3 decay observables is the frequentist i?fit 
framework of the CKMfitter Group [1, 2]; in future work 
we will explore the other analysis schemes outlined by 
the CKMfitter Group [1, 2] and the UTFit Collaboration 
[3-7], such as Bayesian analyses. We review CKMfitter's 
iJfit analysis in sufficient detail to provide sufficient con- 
text for the discussion of our global fit. Complete details 
on the iJfit statistical framework can, of course, be found 
in their original papers [1, 2]. 

Our global fit includes two different types of exper- 
imental observables: (i) results for angular correlation 
coefficients as a function of (binned) electron energy and 
(ii) the neutron lifetime. We consider each of the bin- 
by-bin results for the angular correlation coefficients to 
constitute a separate result. Adopting the notation of the 
iJfit framework, we label each of these experimental mea- 
surements a;exp,i- Each of them are then compared with 
a corresponding theoretical calculation of that quantity, 
2;theo,j- These theoretical calculations arc each a function 
of A^mod model parameters, the set of which we denote 
as {t/mod}- Of these N^^d parameters, Afree < A^mod 
arc experimentally-accessible "free parameters" of the 
model, the set of which we denote as {yfree}- The re- 
maining Acaic = A^mod — Afree "calculated parameters," 
for which there have been no prior experimental mea- 
surements and which are not accessible in the current 



experiments (noting, e.g., second-class currents), must 
be calculated within the context of the model, subject to 
various assumptions. The set of these calculated param- 
eters we denote as {ycaic}- 

The set of experimental observables {xexp} in- 
cludes binned- in-energy measurements of ^exp,i(-Eej) 
and aexp,i(-Be,j) [where we use the subscripts i and j 
to label the particular experiment and the energy bin, 
respectively] and results for the neutron lifetime, Texp.i, 
from different experiments, 

{s^exp} = {^exp,l(£'e,l), ^oxp.l (-£'6,2) j • • • 7 

^exp,l t t^exp,l (i^e,2),..., 

}. (25) 

These are then to be compared, one-by-one, with a cor- 
responding set of theoretical calculations, 

{2;theo(?/mod)} = { Aheo,l (-S'e,!) , ^thco,l (-^6,2) j • • • , 

C*theo,l(-E'e,l), Othco,! (£^e,2),..., 

}, (26) 

which depend on the set of {j/mod} parameters. Under 
the SM, the set of {ymod} parameters would include 

{2/mod} = {A, /2, /s, 92, gs, Vud} , (27) 

though we note that gs does not appear in /3 decay 
observables computed through NLO precision. Conse- 
quently we set (73 = in all that follows; we refer the 
reader to the reviews of Rcfs. [74, 75] for information on 
this quantity. In the next section where we show results 
from example fits for different scenarios, we define for 
each scenario which of the UmodA parameters are to be 
considered a free or calculated parameter. We note a 
subtlety in regards to A: if it is a fit parameter, rather 
than a calculated one, it can be modified by the appear- 
ance of a right-handed coupling emergent from physics 
BSM [33]. Thus the determined value of A from a fit of 
neutron beta decay observables need not be equivalent 
to the calculated value of A = gA/gv- However, present 
lattice calculations of arc of rather poorer precision 
than empirical determinations. The above sets could, of 
course, be trivially expanded to accommodate measure- 
ments of other observables; for example, the {.Tcxp} set 
could include measurements of the neutrino asymmetry 
B and (direct) measurements of the Fierz interference 
term 6bsm (via spectral shape measurements), and/or 
the {ymod} set could include 5bsm- 

As per the usual prescription [22], we define our 
function in terms of a likelihood function £(ymod) for the 
{z/mod} parameter set as 

x'(ymod) = -21n£(y„,od)- (28) 

Following the i?fit framework [1, 2], we define C{ymod) to 
be the product of an "experimental likelihood" function, 
£exp, and a "theoretical likelihood" function, £caic, 

■^(Z/mod) — '£cxp({^cxp}7 {^theo(ymod)}) 

X£calc({2/calc})- (29) 
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The experimental likelihood function is defined in the 
usual way to be the product of the likelihood functions 
for each of the A^exp experimental results, 



■^exp({*^exp}5 {*^theo}) J_ J_ ^exp,i 



(30) 



where, in the ideal case, the individual likelihood func- 
tions are taken to be Gaussian, 



-'exp,i 



: exp 



exp,z 



exp,'; 



(31) 



where cToxp.i denotes the statistical uncertainty of the zth 
experimental result. Of course, the individiial experi- 
mental likelihood functions must account for systematic 
errors, and the formaUsm for the inclusion of such within 
the context of the i?fit framework is described in detail in 
Refs. [1, 2]. However, a detailed discussion of the impact 
of experimental systematic errors is beyond the scope of 
this paper, as the focus of our first paper is on the sta- 
tistical impact of a global fit and the limitations on such 
from theoretical uncertainties. 

In principle, a theoretical likelihood function could 
similarly be defined as the product of likelihood func- 
tions for each of the j/caic.i calculated parameters. 



rcalc({?/calc}) = H -^-Ic,, 



(32) 



i=l 



and under the assumption that the r/caici values are 
Gaussian distributed, the theoretical likelihood would, 
by definition, contribute to the x^- Such a formulation 
might not be appropriate for the treatment of the j/caic,i 
parameters, for which the underlying probability distri- 
butions are certainly not known. However, it may well 
be possible to bound the value of each t/caic,i parame- 
ter on theoretical grounds, such that the parameter may 
rciasonably assume any value over an allowed range of 
[Vcaici - Syca\c,i,yca\c,i + ^J/caici]- It would be highly un- 
likely that the true value of the parameter would fall 
outside of this range. 

Thus, given the lack of knowledge on the underlying 
distributions of the 2/caic,i parameters, the proposal of the 
i?fit scheme [1, 2] is to redefine the function so that 
the theoretical likelihood does not contribute to the x^, 
while the ycaic.i parameters are permitted to vary freely 
within their pre-defined allowed ranges. In particular, 
the x^ is re-defined to be 



Wexp ^ \ 2 

2 V ^ I •^exp,* '^theo,i 

X — / , 



i=l 



exp,i 



-21n£ealc({2/calc}), (33) 



where 



-21n£calc({ycalc}) 



0, Vj/calci e [t/calc,i ± %alc,i] 

00, else 



Thus, under the i?fit scheme, each of the ycaic.i param- 
eters are bounded, but all possible values of the param- 
eters within their pre-defined ranges are treated equally. 
That is, the value of x^ is scanned over the available 
{yfree} parameter space, while the values of the {j/caic} 
parameters are permitted to vary freely over their pre- 
defined ranges at each point in the {t/free} parameter 
space. Thus, the central challenge of such an analysis 
in Refs. [1, 2] is to define the [jjcaicA i (5ycaic,i] allowed 
ranges carefully because: (i) the fit results for the {yhee} 
parameters can be interpreted as valid only if the "true" 
values for the {ycaic} parameters fall within the allowed 
ranges; and (ii) choosing the allowed ranges to be too 
wide (i.e., too conservative) could mask the discovery of 
new physics. 

After construction of the x^) a global fit can be then be 
pursued under two different types of analyses: (i) deter- 
mining values for the SM parameters; and (ii) assessing 
the validity of the SM. 



B. Determining Standard Model Parameters 

Here the goal is neither to assess the validity of the SM 
nor to search for evidence of new physics. Instead, the 
SM is assumed to be valid, and the global fit is employed, 
optimally, to determine values for all of the {^mod} pa- 
rameters. In this case, the minimum value of x^iUmod), 
computed according to Eq. (28), is obtained by allowing 
all of the A^mod parameters to freely vary. The resulting 
minimum value is denoted x^(2/mod)min- Confidence lev- 
els, ^(ymod), on the values of the parameters obtained at 
X^(2/mod)min are calculated according to 



Vivmod) = Prob(Ax'(ymod),iVdof), 



(34) 



where, as usual, 



Ax^(ymod) 



X^(j/mod) - X^(ymod)min, (35) 



and Prob(. . . ) denotes the probability for a value of x^ > 
Ax^(ymod) for A^dof degrees of freedom. 

However, in practice, it may not be possible, or fea- 
sible, to determine values for all of the {ymod} param- 
eters. The {j/caic} parameters, notably the second-class 
couplings, /s and g2, and the induced pseudoscalar cou- 
pling (^3, serve as examples. We let {ya} denote the 
Na < A^niod subset of {t/mod} parameters for which the 
goal is to determine confidence levels; the remainder of 
the Nfj^ = A^niod — parameters are denoted {y^}. Note 
that the set of {j/;^} parameters may not be identical to 
the set of {j/caic} parameters. 

Confidence levels on the {pa} parameter set are then 
determined by first computing at each point in the {ya} 
parameter space the minimal value of the x^ function, 
X^iiUa}; {yti})min, obtained by allowing the {y^} param- 
eters to vary. The minimal value of Ax^ at that point in 
the {ya} parameter space is then computed according to 



Ax'({2/a})=x'({ya};fe})r 



X (2/mod)min- (36) 
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The confidence levels are then obtained from 

V{{ya}) = Prob(Ax'({j/J), iVdof). (37) 

C. Assessing the Standard Model 

Under the minimization scheme just described for the 

determination of SM parameters, the SM is assumed to 
be valid by definition. In principle, a test statistic for 
assessing the validity of the SM would be the value of 

X^(?/mod)min obtained when all of the A^mod arc varied, 
where a confidence level on the SM could be defined as 

P{SM) < Prob(x'(t/mod)min, iVdof) • (38) 

In practice, an assessment of the validity of the SM can 
be obtained via a Monte Carlo calculation according to 
the following scheme. Values for the set of {xoxp} experi- 
mental results are sampled in the Monte Carlo from their 
corresponding set of theoretical expressions, {a;theo}, as- 
suming the fitted values for the {ymod} parameter set. 
For each set of {.Toxp} values, a value is computed, as 
before, by allowing the ymod parameters to vary. This is 
repeated, and from this a (normalized) distribution of x^ 
values, p{x^), is constructed. A confidence level for the 
SM is then deduced from this distribution according to 

^(SM)< J p{x')dx^. (39) 

X^>X^(!/mod)min 

V. EXAMPLE FIT SCENARIOS 

We now illustrate, via several examples of nFitter fits 
to Monte Carlo pseudodata, the impact that simultane- 
ous fits to the energy dependence of the a and A angu- 
lar correlation coeflBcients have on an assessment of the 
validity of the SM and the extent to which theoretical 
uncertainties can limit such an assessment. 



A. Monte Carlo Data Sets and Statistics 

We generated Monte Carlo pseudodata by sampling 

the relevant recoil-order diS'erential distributions for 
measurements of a and A, Eqs. (7) and (15), respec- 
tively, employing the complete expression for h(x) in 
Ref. [56] in the latter case — only the terms of Eq. (17) 
are appreciably nonzero. We generated two different data 
sets. The first data set, which we term our "Standard 
Model" data set, consists of 5 x 10^ simulated events for 
a measurement of o, and a separate data set consisting of 
5 X 10^ simulated events for a measurement of A. Both 
data sets employ the current Particle Data Group av- 
erage value for A = 1.2701 [22] and the CVC value for 
/2 = (np - K„)/2 = 1.8529450 [22], and assume all of the 
other small terms are zero, /a = <?2 = ffs = 0. 



The second data set, which we term our "New Physics" 
data set, again consists of 5 x lO*' simulated events for a 
measurement of a, and 5 x 10^ events for a measurement 
of A. This data set is identical to the earlier one, save 
for the inclusion of a nonzero value for a tensor coupling, 
namely, gr^T = 1-0 x 10~^, close to the strongest em- 
pirical limit on this quantity, which comes from a Dalitz 
study of radiative pion decay [18, 76]. In specific, we 
note the extracted 90% CL limit on Re(eT) [18, 76] can 
be combined with gx = 1.05(35) [33], or gx < 1.4, both 
in the MS scheme at a renormalization scale of 2 GeV, 
to yield -1.5 x 10"^ < grRcieT) < 1.9 x lO'^. The 
most stringent limit on (7sRe(es), which comes from the 
analysis of 0+ 0+ nuclear decays [25], is of a com- 
parable magnitude. However, as can be seen from Eq. 
(13), for approximately equal values of gses and gx^T, 
&BSM is significantly more sensitive to tensor couplings; 
hence, we have decided to illustrate our methods using a 
non-zero tensor coupling exclusively. A summary of our 
input parameters and the resulting values for S, ag, Aq, 
bssM, and t are given in Table I. Note that we calcu- 
late T as per Eq. (20), assuming the central value for the 
superallowcd 0+ ^ 0+ value of T4d = 0.97425(22) [22]. 

Our pseudodata consists of 5 x 10^ events, because 
such is needed for the anticipated level of statistical pre- 
cision in the most ambitious of the next generation of 
decay correlation experiments. Specifically, we would be 
able to determine the value of A in an a measurement to 
0.010% (equivalent to a 0.037% determination of uq) and 
to 0.008% in an A measurement (0.032% determination 
of Aq). The stated goals on a and A in the upcoming 
PERC experiment [20] are to achieve statistical and sys- 
tematic errors on the level of 0.03%. 

The numerical results presented hereafter employ the 
full Eg energy range (i.e., kinetic energies < Tg < Tq). 
Of course, in a real experiment, the lower energy range 
will necessarily be in excess of zero due to hardware 
thresholds and/or analysis cuts. Also, considerations of 
systematics may limit the upper energy range because en- 
ergy loss effects become disproportionately more impor- 
tant as the electron energy increases, see, e.g., Refs. [42- 
44]. Such details would, of course, be included in a global 
fit to actual data; the point of this paper is to illustrate 
the method. 



B. Examples: Fits to the Standard Model Data Set 

For our first example, we consider fits to the Stan- 
dard Model data set. As an illustration of our meth- 
ods, we show our simulated data for a and A as a func- 
tion of Te in Fig. 1, employing 79 10-keV bins from 0- 
790 keV, as the endpoint is Tq = 781.5 keV. We plot 
the "experimental asymmetries" a^^p and A^^p, scaled 
by the nominal electron-energy dependence of the 
asymmetries, where the factor of ^ results from the an- 
gular integral (cos e,^) = | on a hemisphere. These 
experimental asymmetries Oexp and A^^p are calculated 



10 



TABLE I. Summary of parameters for the "Standard Model" and "New Physics" Monte Carlo pseudodata sets. We note under 
CVC that /a = (sp - k„)/2 = 1.8529450 [22]. 



Input Parameters 


A 




/2 


h 


.92 


.93 


3ses 


qt^t 


Standard Model 


PDG: 1.2701 


CVC 



















New Physics 


PDG: 1.2701 


CVC 


{f^p f^n^/2 














1.0 X 10"^ 


Calculated Parameters 






ao 


Ao 


&BSM 


T 






Standard Model 


5.83946 




0.105002 


-0.117495 





885.631 s 






New Physics 


5.83951 




0.104998 


-0.117489 


-0.00522 


885.624 s 








-0.115- 



100 200 300 400 500 600 700 800 
Te [keV] 



-0.11 



-0.115 



-0.12- 



-0.125 



I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 



(b) 



_n 1^1 ■ I ■ I ■ ■ I ■ I ■ I I ■ I ■ I I I ■ ■ I I I ■ ■ ■ I I ■ I ■ I I ■ I ■ I ■ ■ I ■ I 
100 200 300 400 500 600 700 800 
Te [keV] 



FIG. 1. (Color online) Simulated data from the Standard Model data set for acxp/(|/3) [panel (a)] and ^cxp/(|/3) [panel (b)] 
plotted as a function of T^. The solid red line is the result of a simultaneous fit to the a and A data. 



from the simulated data in a manner similar to how ac- 
tual experimental data would be analyzed in a typical 
"forward/backward" asymmetry measurement (see, e.g., 
[77]), where 



7V(cos6'ei. > 0) - A^(cos( 
N {cos 6eu > 0) + A^(cos( 



<0) 
<0) 



1 



ai 



2^1 + &BSM^ 

NicosOe > 0) - 



' 3 

iV(cos 



<0) 



^exp 



N{cos9e > 0) +N{cosl 
2^1 



<0) 



&BSM- 



(40) 



(41) 



Sensitivity to 6bsm from A{E) and a{E) have been pre- 
viously considered by Refs. [33, 39, 79]. In a real experi- 
ment the effects of 0{a) radiative corrections [78] would 
have to be removed to interpret Aexp in terms of the sim- 
ple theoretical expressions we employ, noting Eqs. (8) and 
(18), in our fits. We avoid this now for simplicity, and 
we are able to do so because said correction incurs no 
additional hadronic uncertainty. Moreover, for similar 
reasons we drop the 02 term from our fits as well; they 
are simply trivially small. The fits shown in Fig. 1 are 
the result of a simultaneous fit to the a and A data, in 
which {xcxp} — {o.,A}, noting {a, A} is shorthand for 
the complete set of the binned-in-energy results for Oexp 
and ^cxp, and {t/a} = {A}. We fix /2 to its CVC value 



and set all second-class couplings to zero, so that &bsm 
vanishes. As a validation of our methods, the fit result 
for A = 1.27009(8) agrees with the input value to within 
— O.lcr with a Xmin/-^dot — 135.3/157, yielding a perfectly 
acceptable Prob(x^ > Xmin) = O-^^- 

Relaxing the assumption that second class currents 
are zero, we apply the _Rfit scheme to a fit in which 
{Va} — {A}, and /a and §2 comprise the {y^} param- 
eter set, which are then permitted to vary simultane- 
ously over some prescribed range, as per the prescription 
discussed in Sec. IV B. Of the other {?/mod} parame- 
ters, /2 is again fixed to its CVC value, and 6bsm is 
fixed to zero. The resulting 68.3% CL on A for differ- 
ent assumptions on the permitted theory ranges for /a 
and/or 32 are compared in Table II. Note that we deter- 
mine a 68.3% CL as per the requirement Ax^({2/a}) = 
X^iiVa}] {y^Dmin " (2/mod)mm = 1, whcrc in this case 
{Va} = {A} and {y^} = {/3,ff2}- Referring to Table 
II, unless g2 can be constrained to 0(0.1), theory uncer- 
tainties in §2 would limit the precision to which A can 
be extracted from experiments aiming to measure a and 
A to the level of 0.03%. Even at this level, the range 
of the 68.3% CL on A is ~ 50% larger than the case in 
which second class currents are taken to be exactly zero. 
In contrast the fits are almost completely insensitive to 
the value of /a; this is because the latter appears only in 
the e/Rx terms. 

Alternatively, in the absence of a theory bound on g2 , 
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TABLE II. Fitted values for A from a simultaneous fit to the 
{a, A} Standard Model data set under the inclusion of the 
indicated theoretical uncertainties in fs and 32. The fitted 

values for A arc defined by the location of the overall Xmim 
while the 68.3% CL is defined by A^^ = 1. An asymmetric 
CL about the Xmin value is indicated by asymmetric upper 
(+) and lower (— ) error bars. Recall that the input value for 
A = 1.2701 [22]. 



/a Range 


P2 Range 


Fitted A 








1.27009(8) 




[-0.1,0.1] 


[-0.1,0.1] 



1.27036(+7)(-17) 
1.27009(8) 




[-0.5,0.5] 


[-0.5,0.5] 



1.27056(31) 
1.27007(+ll)(-7) 


[-0.1,0.1] 
[-0.5,0.5] 


[-0.1,0.1] 
[-0.5,0.5] 


1.27035(+9)(-17) 
1.27054(+32)(-31) 



TABLE III. Fitted values for A and 32 from two-parameter 
simultaneous, {a, A}, and individual, {a} or {A}, fits to the 
Standard Model data set. The fitted values for A and 52 are 
defined by the location of the overall Xmim while the 68.3% 
CL is defined by Ax^ = 2.30. We use the inputs A = 1.2701 
[22] and (72 = 0. The column "Prob" indicates the Prob(x^ > 

Xmin)- 



Fit 


Fitted A 


Pull 


92 


Xmin/^dof 


Prob 


{a, A} 
{a} 

{A} 


1.27056(48) 
1.294(40) 

1.2708(14) 


0.97(7 
0.60(T 

O.SOcr 


0.174(171) 
6.7(11.3) 

0.28(65) 


133.0/156 
54.4/77 
77.8/77 


0.91 
0.98 
0.45 


Fit 


Correlation Coefficient PA92 






{a, A} 

{a} 

{A} 


0.97 
1.00 
0.99 











one could fit directly for A and 52- In what follows we 
neglect the small fs contribution, setting = 0, which 
is reasonable, as we show in Table IL Results from a 
fit in which {ya} = {A, (72}, with no parameters, are 
shown in Table III, where the errors on A and 52 are de- 
fined by Ax^ 2.30, i.e., the 68.3% CL for a joint fit 
of 2 free parameters [22]. As one would expect, the er- 
ror on A from a two-parameter simultaneous fit to the 
{a, A} data set is a factor of ~ 6 larger than that from 
a single-parameter fit for A alone; and, it is worthy to 
note that 92 can be determined from such a fit to 0(0.2). 
However, what is more interesting are the errors on A 
and (72 extracted from a fit to the {a} data set alone: the 
error on A is a factor of ~ 300 larger than that from a 
single-parameter fit for A alone, and the error on 52 is of 
O{10). The origin of this effect is clear: g2 appears in oi 
only in expansion for a, Eq. (8), namely via the combi- 
nation 2(A — A^)(l +52)- Therefore, A and 52 are directly 
correlated in a small recoil-order term in a, with only the 
fitted value of ao ultimately limiting the Ax^ range and, 
hence, the A and 52 errors. In Table III, we also show 
the "Pull," which we define as Pull = {xf^t — a;input)/(o'fit) 
for a parameter x. In the event of asymmetric errors, we 



TABLE IV. Fitted values for A from a simultaneous fit to the 
{a. A} Standard Model data set under different assumptions 
on the range of CVC breaking for the value of /2 . The fitted 

values for A are defined by the location of the overall Xmin; 
while the 68.3%i CL is dc^fincd by Ax^ = 1. Wo usc^ the input 
A = 1.2701 ;22]. 



/2 Range 


Fitted A 


CVC Exact 


1.27009(8) 


±1% 


1.27009(9) 


±2% 


1.27011(11) 


±5% 


1.27016(15) 



average the two errors to form aat- For completeness, the 
correlation coefficients p\g^ from the fits are also given. 
Thus, this example succinctly illustrates the necessity of 
a stringent theoretical bound on the value of (?2 for an 
interpretation of a measurements in the context of an 
assessment of the V — A structure of the SM. 

As a final example of fits to our Standard Model data 
set, we consider the implications of CVC breaking on the 
value of the weak magnetic form factor /2, by defining 
{Ua} = {A} and {y^} = {/2} under different assump- 
tions on the permitted theory range for /2. The results 
are summarized in Table IV. As can be seen there, the 
impact of a ±2% breaking on /2 is comparable to an 
0(0.1) uncertainty in 52- 



C. Examples: Fits to the New Physics Data Set 

As our second example, we consider fits to the New 
Physics data set. As can be seen in Table I, the values of 
ao and Aq in the New Physics data set differ from their 
values in the Standard Model data set by only 0.004% 
and 0.005%, respectively. Therefore, the impact of any 
new physics from scalar and tensor interactions on the 
measured values of a and A will be via a "dilution" to 
the experimental asymmetries a^^p and A^^p from a non- 
zero Fierz term 6bsm appearing in the denominators of 
Eqs. (40) and (41), respectively. Accordingly, we now 
expand our {j/mod} parameter set to include &bsm- 

Results from a single parameter fit in which {ya} = {A} 
only and an empty {y^j,} parameter set are shown in Table 
V. As can be seen, these fits still yield excellent values for 
xj^jii^; however, not siirprisingly, there arc significant pulls 
on the fitted values for A from the input value. The origin 
of these pulls can be easily understood (and, indeed, was 
first noted by Ref. [33]) by inspecting the functional forms 
for Oexp and A^xp in terms of a and A. For small 6bsM) 
CKexp ~ Q;(l — bBSMTrie/Ee) with a = a or A, whereas the 
slope of the e/Rx oc nie/E,, dependence of ai and A is 
negative in A and nearly linear; the rrie/Ee contributions 
to ttexp and ^exp from 6bsm and A are of the same sign 
and slope. Therefore, the presence of a non-zero 6bsm 
would not result in poor values in fits to the a and A 
energy dependence. 



TABLE V. Fitted values for A from simultaneous, {a, A}, and 
individual, {o} or {A}, fits to the New Physics data set, in 
which A was the only free parameter. The fitted values for A 

are defined by the location of the overall Xmin in each case, 
while the 68.3% CL is defined by Ax^ = 1. We use the 
input A = 1.2701 [22]. The column "Prob" indicates the 

Pr0b(x'' > Xmin)- 



Fit 


Fitted A 


Pull 


Xmin/-^dof 


Prob 


{a, A} 


1.27115(8) 


13.1(7 


153.6/157 


0.56 


{a} 


1.27135(13) 


9.6a 


75.0/78 


0.57 


{A} 


1.27103(10) 


9.2(7 


74.7/78 


0.59 



TABLE VI. Fitted values for A from simultaneous, {a, A}, 
and individual, {a} or {A}, fits to the New Physics data 
set, under the inclusion of the indicated theoretical uncer- 
tainties in /a and g2. The fitted values for A are defined 

by the location of Xmim while the 68.3% CL is defined by 
Ax^ = 1. An asymmetric CL about the Xmin value is indi- 
cated by asymmetric upper (+) and lower (— ) error bars. We 
use the input A = 1.2701 [22]. The column "Prob" indicates 

the Prob(x^ > Xmin)- 



Fit 


fa € [-0.1,0.1] and 92 € [- 
Fitted A Pull 


-0.1,0.1] 

Xmin /-^dof 


Prob 


{a, A} 

{a} 

{A} 


1.27088(+19)(-7) 6.0(7 
1.27170(+13)(-83) 3.3(7 
1.27124(+10)(-27) 6.2(7 


151.9/157 
75.0/78 
73.9/78 


0.60 
0.58 
0.61 


Fit 


fs € [-0.5, 0.5] and 32 € [- 
Fitted A Pull 


-0.5,0.5] 

Xmin /-^dof 


Prob 


{a, A} 

{a} 

{A} 


1.27072(+30)(-32) 2.0.7 
1.27311(+15)(-365) ~0 
1.27209(+13)(-52) 3.2a 


151.4/157 
74.8/78 
71.8/78 


0.61 
0.58 
0.68 



Next, we again consider the implications of theoretical 
uncertainties in /a and 32 as relevant for an extraction of 
A in the presence of new physics. The {ua} parameter set 
still includes A only, while now the {y^} parameter set 
includes /a and 52, which are then permitted to vary over 
a particular range. The resulting 68.3% confidence levels 
on A extracted from a simultaneous fit to the {a, A}, as 
well as to the individual {a} or {A}, New Physics data 
set, for different theory ranges on /a and 92 are compared 
in Table VI. It is interesting to note that the effects of 
new physics in an a measurement considered in isolation 
could potentially be obscured by second class currents, 
whereas there is less of an effect for an A measurement 
considered in isolation or in a combined fit to a and A 
data. This effect derives from the manner in which g2 
appears in a. 

Results from the two-parameter fits, with {Ua} = 
{A, 6bsm} and an empty {y^} parameter set, are sum- 
marized in Table VII. It is worth noting that a com- 
bined {a, A} fit at this level of precision has the po- 
tential to constrain 6bsm at 68.3% CL to the level of 
3 X 10~^, although such a fit would also offer sig- 
nificantly less (by a factor of ^ 8) less sensitivity to 
A. For completeness, we also list the correlation coef- 
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TABLE VII. Fitted values for A and 6bsm from simultaneous, 
{a. A}, and individual, {a} or {^4}, fits to the New Physics 
data set. The fitted values for A and 6bsm are defined by 

the location of the overall Xmin; while the 68.3% CL is de- 
fined by Ax^ = 2.30. We use the inputs A = 1.2701 [22] 
and &BSM = —0.00522. The column "Prob" indicates the 

Prob(x^ > Xmin)- 



Fit 


Fitted A 


Fitted 6bsm 


Xmin/-^dof 


Prob 


{a, A} 


1.27011(65) 


-0.0051(31) 


147.5/156 


0.67 


{a} 


1.27052(113) 


-0.0037(50) 


73.8/77 


0.58 


{A} 


1.27014(86) 


-0.0045(44) 


72.2/77 


0.63 


Fit 


Correlation Coefficient pxb 






{a, A} 


0.98 








{a} 


0.99 








{A} 


0.98 









0| 1 1 1 1 1 1 1 1 1 1 j 1 1 1 1 1 1 1 r 




1.271 



FIG. 2. Impact of uncertainties in /s and g2 on the allowed 
(A, feBSin) parameter space from a two-parameter simultane- 
ous fit to the {a, A} New Physics data set. The bands indi- 
cate the 68.3% CL allowed regions defined by Ax^ = 2.30 for a 
joint fit of 2 free parameters, for the indicated uncertainties in 
/a and 52. The input values A = 1.2701 and 6bsm = -0.00522 
are indicated by the dashed lines. 



ficients pxi, for these fits. Finally, we perform fits with 
{Va} = {A,6bsm} and {y^} = 1/3,92}, with /s and 52 
permitted to vary over different particular ranges. The 
results of these fits, which demonstrate the impact that 
uncertainties in /s and 52 have on the allowed (A, 6bsm) 
parameter space, are shown in Fig. 2. Indeed, as one can 
see, the allowed (A, 6bsm) parameter space is broadened 
significantly by the inclusion of theoretical uncertainties 
in /a and 32- As per our analysis reported in Table II we 
note the uncertainty in g2 is of greater impact. 



D. Example: Impact of the Neutron Lifetime to an 
Assessment of the SM 

Finally, we consider the impact that future measure- 
ments of the neutron lifetime to a precision of 0.1 s, con- 
sidered together with measurements of the a and A an- 
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TABLE VIII. Results for Xmin from a simultaneous fit for A to 

the {a, A, t} New Physics data sot, for different assumed pre- 
cisions in the measurement of the neutron hfetime, assuming 
92 = 0. 



Tn Error 


Xmin /-^dof 


Prob(x^ > Xmin) 


1.00 s 


155.1/158 


0.55 


0.50 s 


159.3/158 


0.45 


0.25 s 


174.3/158 


0.17 


0.20 s 


183.8/158 


0.078 


0.15 s 


200.9/158 


0.011 


0.10 s 


233.0/158 


8.2 X 10"'' 


o.or s 


2()3.1/io<S 


2.1) X 10 ' 


TABLE IX. Results for Xmin from a simultaneous fit for A and 


g2 to the {a, A, r 


} Now Physics data set, for 


different assumed 


precisions in the measurement of the neutron lifetime, now for 


92 € [-0.025,01. 


The first set of values ignore the role of 92, 


as per usual procedures, in the theoretical formula for t; the 


second set include it. 




Tn Error 


Xmin/-^''dof 


Prob(x^ > Xmin) 


1.00 


154.3/157 


0.55 


0.50 


158.0/157 


0.46 


0.25 


171.1/157 


0.21 


0.20 


179.5/157 


0.11 


0.15 


194.5/157 


0.023 


0.10 


222.6/157 


4.5 X 10"" 


0.07 


249.1/157 


3.9 X 10"® 


1.00 


154.4/157 


0.54 


0.50 


158.3/157 


0.46 


0.25 


172.1/157 


0.19 


0.20 


180.9/157 


0.093 


0.15 


196.7/157 


0.017 


0.10 


226.3/157 


2.5 X 10"^ 


0.07 


254.3/157 


1.4 X 10"® 



gular correlation coefficients to a precision of ~ 0.03%, 
will have on an assessment of the validity of the V — A 

structure of the SM. We illustrate this within the con- 
text of our New Physics data set. To do so, we expand 
our {a;exp} data set to {a, A, r}, where we take r to be 
a single data point whose value is the calculated value of 
T = 885.624 s for the parameters of the New Physics data 
set, as per Table I, and consider various uncertainties in 

T. 

We then perform a simultaneous fit to the {xoxp} = 
{a, A,r} New Physics data set in which {ya} — A, i.e., 
with only a single free parameter. We use the CVC value 
for /2 and set /s = 52 to zero. As we have already noted, 
a simultaneous fit to the {a. A} data set only yields ex- 
cellent values for Xmin — albeit, with significant pulls 
on the fitted values for A. The neutron lifetime is rela- 
tively insensitive to new scalar or tensor physics, because 
such new physics only enters quadratically in the 5 pa- 
rameter, noting Eq. (20). Therefore, we would expect a 
simultaneous fit to an {a. A, r} data set in the presence 
of new scalar and tensor interactions to return a poor 



TABLE X. Results for Xmin from a simultaneous fit for A and 

52 to the {a, A, r} New Physics data set, for different assumed 
precisions in the measurement of the neutron lifetime, now for 
92 G [—0.1,0.1]. The first set of values ignore the role of 52, 
as per usual procedures, in the theoretical formula for t; the 
second set include it. 



Tn Error 


Xmin /-^dof 


Prob(x^ > Xmin) 


1.00 


152.7/157 


0.58 


0.50 


155.0/157 


0.53 




1 CQ Q /I 




0.20 


168.6/157 


0.25 


0.15 


178.0/157 


0.12 


0.10 


195.7/157 


0.019 


0.07 


212.5/157 


0.0021 


1.00 


152.9/157 


0.58 


0.50 


155.9/157 


0.51 


0.25 


166.5/157 


0.29 


0.20 


173.3/157 


0.18 


0.15 


185.4/157 


0.060 


0.10 


208.1/157 


3.9 X 10"^ 


0.07 


229.6/157 


1.4 X 10"" 


TABLE XI. Results for Xmin from a simultaneous fit for A and 


92 to the {a, A, t} New Physics data set, for different assumed 


precisions in the measurement of the neutron lifetime, now for 


92 G [-0.5,0.5]. 


As in Table IX the first set of values ignore 


the role of 92 in 


the theoretical formula for t; the second set 


include it. 






Tn Error 


Xmin/-^dof 


¥roh{x^ > Xmin) 


1.00 


152.1/157 


0.60 


0.50 


153.0/157 


0.58 


0.25 


154.3/157 


0.55 


0.20 


154.7/157 


0.54 


0.15 


155.0/157 


0.53 


0.10 


155.3/157 


0.52 


0.07 


155.4/157 


0.52 


1.00 


152.4/157 


0.59 


0.50 


154.3/157 


0.55 


0.25 


158.6/157 


0.45 


0.20 


160.3/157 


0.41 


0.15 


162.2/157 


0.37 


0.10 


164.3/157 


0.33 


0.07 


165.6/157 


0.30 



value for Xmin- Indeed, the results of such an analysis 
are shown in Table Vlll, where we show the fitted val- 
ues for A, the Xmin from the fit, and the probability for 
that x^jjj value, for different assumed experimental er- 
rors on the lifetime. As can be seen, the probability for 
the validity of the SM, which we quantify via the statistic 
P{SM) < Prob(x^ > Xmin). as per Section IV C, in the 
presence of new tensor physics at the level of grer ~ 10"^ 
would become < 10"^ if the neutron lifetime were mea- 
sured to ~ 0.1 s or better in concert with ~ 0.03% mea- 
surements of a and A. Under this scenario, a precision in 
the neutron lifetime of 0.07 s would yield a probability of 
~ 3 X 10"^, which is slightly more stringent than the re- 
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quircmcnt for a 5(t result, i.e., a probability of 5.7 x 10~ . 

Finally we turn to an examination of the impact of a 
non-zero second-class coupling 92 on the ability to fal- 
sify the V — A law of the SM. It has been the usual 
procedure to ignore certain recoil corrections in the de- 
termination of T, as per Eq. (20), but we expect that 
a non-zero value of g2 could be important in this con- 
text, so that we include the recoil correction in 52, as 
per Eq. (24), as well. We perform a simultaneous fit to 
the {Xexp} = {(i,A,t} New Physics data set in which 
{Va} = ^ and {y^} = {g2}, where 52 is permitted to 
vary over different particular ranges, with /2 equal to its 
CVC value and /s = 0. The empirical limits on 52 are 
markedly weaker than the existing direct theoretical esti- 
mate, noting 52 = -0.0193 ±0.0067 [73] with A = 1.2701 
[22] , so that we perform simultaneous fits using g2 in the 
following ranges: 92 € [-0.025,0], 32 G [-0.1,0.1], and 
92 S [—0.5, 0.5]. The fit results, as well as the determined 
abilities to falsify the SM, as a function of the error in the 
determined neutron lifetime, arc shown for these ranges 
of 52 in Table IX, Table X, and Table XI, respectively. 
Nonzero values of 52 impact the ability to falsify the SM 
in every case, and the inclusion of the recoil corrections 
to T are also of importance. In the last case, in which 
92 € [—0.5, 0.5], the ability to falsify the SM with improv- 
ing precision in the neutron lifetime has been completely 
eroded. Evidently it is important to determine 92 to the 
greatest accuracy possible in order to be able to falsify 
the V - A law of the SM. 



VI. SUMMARY AND CONCLUSIONS 

In summary, we have developed a maximum likelihood 
statistical framework, which we term nFitter, in which 
we make simultaneous fits to various neutron ^ decay 
observables. Although a number of global fits to /3 decay 
data have previously been developed [12, 35], the novel 
approach embedded in our technique is that simultaneous 
fits to the energy dependence of the angular correlation 
coefficients allow for a robust test of the validity of the 

V — A structure of the SM in the presence of theoretical 
uncertainties, whereas fits based on integral quantities 
do not. To our knowledge ours is the first study of the 
quantitative ability to falsify the SM and particularly the 

V — A law, after the manner of Refs. [1-7], in the context 
of neutron beta decay observables. 

Our study has consisted of fits to the a. and A angu- 
lar correlation coefficients, as well as to the value of the 
neutron lifetime. We believe that studies of the B an- 
gular correlation coefficient, as well as of 6bsm through 
the electron energy spectrum in /3 decay, will offer impor- 
tant complementary information; and such, as well as any 
additional, concomitant theoretical uncertainties, can be 
incorporated in our analysis framework as well. In our 



current study we have focussed on the role of second-class 
current contributions, most notably on the impact of a 
non-zero g2 coupling, on the ability to identify physics 
BSM. In the course of developing our analysis procedure, 
we have discovered that certain recoil effects to the neu- 
tron lifetime, contrary to the usual view [48], can have 
an impact on our fit results. Moreover, the precise form 
of the recoil corrections depends on experimental details, 
revealing that the corrections change with a finite experi- 
mental acceptance, such as in experiments which extract 
a value of the lifetime from measurements of the decay 
electrons and/or protons. Such considerations warrant 
further detailed study. 

We have explicitly shown that it is possible to dis- 
cover physics BSM, at 5(t significance, in neutron /? de- 
cay observables using experiments which are currently 
planned or under construction. This is subject to the 
following conditions; namely, that (i) tensor interac- 
tions are not much smaller than the constraint which 
emerges from the Dalitz analysis of pion radiative (3 de- 
cay [18, 76], (ii) the value of 92 can be sharply restricted, 
and (iii) results of 0.03% precision can be realized for a 
and A, in concert with a sub-O.ls determination of the 
neutron lifetime. In our study we have assumed that 
grReier) = 0.001, so that using gr = 1.05(35) [33] and 
noting Re(eT) ^ '^'^/■^bsM' ^^^^ would be commensurate 
with the appearance of physics BSM at an energy scale 
of at least Absm ^ 5TeV [31]. We note that existing di- 
rect limits on tensor couplings from nuclear /3 decay are 
much weaker than those from radiative pion decay [18]; 
perhaps new physics effects could be different in pion 
and neutron decays. We note, however, that under the 
assumption that BSM effects appear at energies in excess 
of Absm such effects can only occur from operators be- 
yond mass-dimension six and ought be suppressed. Our 
current analysis framework is also suitable to the discov- 
ery of new scalar interactions as well. 

We have shown that theoretical uncertainties in g2 can 
mitigate the gains made in falsifying the SM through the 
inclusion of precision r results. We thus advocate for 
a determination of g2 using lattice gauge theory tech- 
niques; we suppose that lattice measurements of /2 and 
/s in neutron decay would be useful, too. These con- 
siderations are quite independent of how information on 
bssM is determined. 
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